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Within hydrodynamic approach, we study the Bose-Einstein correlation of identical pions by taking 
into account both event-by-event fluctuating initial conditions and continuous pion emission during the 
whole development of the hot and dense matter formed in high-energy collisions. Important deviations 
occur, compared to the usual hydro calculations with smooth initial conditions and a sudden freeze- 
out on a well defined hypersurface. Comparison with data at RHIC shows that, despite rather rough 
approximation we used here, this description can give account of the bit dependence of Rl and R s and 
produces a significant improvement for R with respect to the usual version. 



PACS numbers: 25.75.-q, 24.10.Nz, 24.60.-k 

Introduction - When describing ultra-relativistic heavy- 
ion collisions in hydrodynamic approach^, a simple pic- 
ture has been extensively adopted. It is usually consid- 
ered that, after the initial very complicated interaction 
between two incident nuclei, at some early instant of time 
a local thermal equilibrium is attained. Such a state is 
usually described in terms of a set of highly symmetric 
and smooth distributions of velocity and thermo- 
dynamical quantities. These are the initial conditions 
(IC) for the hydrodynamic equations, which must be 
complemented with reasonable equations of state (EoS). 
Then, as the thermalized matter expands, the system 
gradually cools down and, when the temperature reaches 
a certain freeze-out value T/ , it suddenly decouples. 
Every observed quantity is then computed on the hyper- 
surface T — Tf a . For instance, the momentum distribu- 
tion of the produced hadrons are obtained by using the 
Cooper-Frye integral|2( extended over this hypersurface. 

Though operationally simple, and actually useful for 
obtaining a nice comprehension of several aspects of the 
phenomena, such a scenery is clearly highly idealized 
when applied to finite- volume and finite-lifetime systems 
as those formed in high-energy heavy-ion collisions. In 
this letter, we examine modification of two ingredients of 
such a description, namely, i ) effects of fluctuations in 
the IC; and ii) consequences of continuous emission 
(CE) of particles, regarding two-pion correlation. 

The identical-particle correlation, also known as HBT 
effect HQ is a powerful tool for probing geometrical sizes 
of the space-time region from which they were emitted. 
If the source is static like a star, it is directly related 
to the spatial dimensions of the particle emission source. 
When applied to a dynamical source, however, several 
non-trivial effects appear 0,0], reflecting its time evolu- 
tion as happens in high-energy heavy-ion collisions. Be- 
ing so, the inclusion of IC fluctuations and of the contin- 
uous emission may affect considerably the so-called HBT 
radii, because both of them modify in an essential way 
the particle emission zone in the space-time. 

The usual symmetric, smooth IC may be understood 



as corresponding to the mean distributions of hydrody- 
namic variables averaged over several events. However, 
since our systems are not large enough, large fluctuations 
varying from event to event are expected. In previous 
publications 0, @ , we showed that indeed the effects of 
these fluctuations on the observed quantities are sizeable 
and moreover in||, we compared the rapidity distribu- 
tions of pions and p — p and showed that the average 
multiplicicity of pions decreases if we consider the event- 
by-event fluctuations, in comparison with the one given 
by the averaged ic[t|. Concerning two-pion correlation, 
as the IC in the event-by-event base often show small 
high-density spots in the energy distribution, our expec- 
tation is that these spots manifest themselves at the end 
when particles are emitted, giving smaller HBT radii. 

As for the decoupling process, it has been proposed [Tol 
llH an alternative picture where the emission occurs not 
only from the sharply defined freeze-out hypersurface, 
but continuously from the whole expanding volume of the 
system at different temperatures and different times. Ac- 
cording to this picture, the large-transverse-momentum 
(fcx) particles are mainly emitted at early times when 
the fluid is hot and mostly from its surface, whereas the 
small- kx components are emitted later when the fluid is 
cooler and from larger spatial domain. Mostly by using a 
simple scaling solution, we showed in the previous papers 
that this picture gives several nice results, namely, i) CE 
enhances the large-m^ component of the heavy-particle 
(p, A, H,fi, ...) spectra, ii) it gives a concave shape for 
the pion to^ spectrum even without transverse expan- 
sion of the fluid[H[l3, Hi) it can lead to the correct hy- 
peron production ratios and spectrum shapes with con- 
ceptually reasonable choice of parameters |lflll2lll3| . and 
iv) it reproduces the observed mass dependence of the 
slope parameter T Concerning HBT correlation, we 
showed 15], within the same approximation, that whereas 
the so-called side radius is independent of the average 
&t 7 the out radius decreases with < kx > , because of 
the reason mentioned above. This behavior is expected 
to essentially remain in the general case we are going to 



2 



discuss below and shown by data|l6j. 
Initial Conditions - In order to produce event-by-event 
fluctuating IC, we use the NeXus event generator fisj. 
based on Gribov-Regge model. Given the incident nu- 
clei and the incident energy it produces the energy- 
momentum tensor distribution at the time \Jt 2 — z 2 — 
1 fm, in event-by-event basis. This, together with the 
baryon- number density distributions, constitute our fluc- 
tuating IC. The strangeness has not been introduced in 
the present calculations. As mentioned in the Introduc- 
tion, we understand that the usual symmetric, smooth 
IC may be obtained from these by averaging over many 
events. In Fig. 1, we show an example of such an event 
for central Au+Au collision at 130 A GeV, compared 
with an average over 30 events. As can be seen, the 
energy-density distribution for a single event (left), at 
the mid-rapidity plane, presents several blobs of high- 
density matter, whereas in the averaged IC (right) the 
distribution is smoothed out, even though the number of 
events is only 30. In the results below, we show that the 
fingerprints of such high-density spots remain until the 
frezeout stage of the fluid, giving smaller HBT radii. 




x (fm) x (fm) 

FIG. 1: Examples of initial conditions for central Au+Au 
collisions given by NeXus at mid-rapidity plane. Energy den- 
sity is plotted in units of GeV/fm 3 . Left: one random event. 
Right: average over 30 random events (corresponding to the 
smooth initial conditions in the usual hydro approach). 

Hydrodynamic Equations - The resolution of the hydro- 
dynamic equations deserves a special care, since our ini- 
tial conditions do not have, in general, any symmetry 
nor they are smooth. We adopt the recently developed 
SPheRIO code|3|, based on the so-called smoothed- 
particle hydrodynamics (SPH), first used in astrophysics 
and which we have adapted for nuclear collisions 19], a 
method flexible enough, giving a desired precision. The 
peculiarity of SPH is the use of discrete Lagrangian co- 
ordinates attached to small volumes ("particles") with 
some conserved quantities. Here, we take the entropy 
and the baryon number as such quantities. Then, the 
entropy density, for example, is parametrized as 

s(x, t) = -^Ui VF(x - Xi{t); h) , (1) 

i 

where W(x — Xj(t); h) is a normalized kernel; Xj(t) is the 
i-th particle position, so the velocity is v j = dxj j dt ; h is 
the smoothing scale parameter, and we have 



S = Jd 3 x 7*(x,t) 
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Observe that, since an entropy Vi is attached to the i-th 
SPH particle, and so is a baryon number, the total en- 
tropy S and the baryon number N are automatically con- 
served. By rewriting the usual energy-momentum con- 
servation equation du T^ — 0, we get a set of coupled 
ordinary equations 19J 
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V i W(x i -x i ;/i). (3) 



Equations of State - For equations of state, we consider 
ones with a first-order phase transition between QGP 
and a hadronic resonance gas, with baryon number con- 
servation taken into account. In the QGP, we consider 
an ideal gas of massless quarks (it, d, s) and gluons, with 
the bag pressure B taken to give the critical temperature 
T c = 160 MeV at zero chemical potential. In the reso- 
nance gas phase, we include all the resonances below 2.5 
GeV, with the excluded volume effect taken into account. 

Two-Pion Correlation - We assume that all pions are 
emitted from a chaotic source and neglect the resonance 
decays. It is argued [2^ that, since resonance decays 
contribute to the correlations with very small q values 
(q ^ q-min , where q m %n is the minimum measureable q) , 
the experimentally determined HBT radii are essentially 
due to the direct pions. Then the correlation function is 
expressed in terms of the distribution function f(x, k) as 
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(4) 



7(0,fc 1 )/(0,fc 2 ) 

where K = (k\ + fca)/2 and q = {k\ — k%) and ki is the 
momentum of the ith pion. Usually 

I(q,K) = (a+a k2 )= [ da^K" f(x,K)e l " x . (5) 

JTjo 

In SPH representation, we write I(q, K) as 

^^ = Err^ e ^'/(^), (6) 

where the sumation is over all the SPH particles. In the 
Cooper-Frye freezeout, these particles should be taken 
where they cross the hyper-surface T = X/ and rij^ is 
the normal to this hyper-surface. 

Continuous Emission Model - In CEM0|i it is assumed 
that, at each space-time point x^, each particle has a 
certain escaping probability 



V{x, k) = exp 



pavdr' 



(7) 
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due to the finite dimensions and lifetime of the thermal- 
ized matter. The integral above is evaluated in the proper 
frame of the particle. Then, the distribution function 
f(x, k) of the expanding system has two components, one 
representing the portion of the fluid already free and an- 
other corresponding to the part still interacting, i.e., 
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f(x,k) = ffree(x,k) + f int (x,k) . 

We may write the free portion as 

ffree(x, fc) = Vf(x, k) . 

The integral JSJ is then rewritten in CEM as 
I(q,K) = [ da ll K»f free (x ,K)e i *' 

J CO 

+ J a*xd li [K»f /ras (x,K)] 



(8) 
(9) 
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where the surface term corresponds to particles already 
free at the initial time. 

The problem of this description is its complexity in 
handling because, V depends on the momentum of the 
escaping particle and, moreover, on the future of the 
fluid as seen in eq.JJJ. In order to make the compu- 
tation practicable, here we first take V on the average, 
i.e., V{x,k) =>< V(x,k) >= V(x) . Then, approximate 
linearly the density p(x') = as(x') in eq.Q. Thus, 



V{x, k) V{x) = exp 



(11) 



\ds/dr\ 
where k = 0.5 a < av >. 

Now, ea. (|10|l is translated into SPH language, by com- 
puting the sum © not over T = Tf a but picking out SPH 
particles according to this probability, with rij^ pointing 
to the 4-gradient of V. Thus, our approximation includes 
also emission of particles of any momentum, once a SPH 
particle has been chosen. However, since our procedure 
favors emission from fast outgoing SPH particles, because 
p decreases faster and so does s in this case making V 
larger, we believe the main feature of CEM is preserved. 

Results - We first assume sudden freezeout (FO) at 
Tf = 128 MeV. This temperature was previously found 
by stud ying the energy dependence of kaon slope param- 



eter T*|21|. It has been also shown[8j that T* is not 
sensitive to IC fluctuations. 

In Fig[21 we compare Ci averaged over 15 fluctuat- 
ing events with those computed from the averaged IC 
(so, without fluctuations). One can see that the IC fluc- 
tuations are reflected in large fluctuations also in the 
HBT correlations. When averaged, the resulting C2 are 
broader than those computed with averaged IC, so giving 
smaller radii as expected. Also the shape of C2 changes. 
We plot the tot dependence of HBT radii, with Gaussian 
fit of C* 2 , in FigO| together with RHIC data[H E3 and 
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FIG. 2: Correlation functions from fluctuating IC and aver- 
aged IC. Sudden freeze-out is used here. The rapidity range 
is —0.5 < Y < 0.5 and q , s ,i which do not appear in the 
horizontal axis are integrated over < q , s ,i < 35MeV. 

results with CEM. It is seen that the smooth IC with sud- 
den FO makes the tut dependence of R flat or even in- 
creasing, in agreement with other hydro calculations [23) 
but in conflict with the data. The fluctuating IC make 
the radii smaller, especially R , without changing the 
TOT-dependence. 

Let us now consider CEM. In this paper, we estimated 
k as being 0.3, corresponding to < av >— 2fm 2 . In 
FigfU we show the charged tot distribution to ensure 
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FIG. 3: HBT radii and the ratio R /R s for sudden freeze-out 
(FO) and CE. 1 stands for averaged IC and 2 fluctuating IC. 
Data are from flH T7|: (tt+ + tt")/2. 
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FIG. 4: charged-particle pr distributions in CEM for the most 
central Au+Au at 130A GeV. The data are from STARp^. 

that the estimate above does reproduce correctly these 
data. Now, look at the mr dependence of the HBT radii, 
shown in FigOU Comparing the averaged IC case, with 
CEM (CE1), with the corresponding frezeout (FOl), one 
sees that, while Rl remains essentially the same, R s de- 
creases faster and as for R Q , it decreases now inverting 
its rriT behavior. The account of the fluctuating IC in 
addition (CE2) makes all the radii smaller as in FO case, 
obtaining a nice agreement with data for Rl and R s , and 
improving considerably the results for R a with respect to 
the usual hydro description. 

Conclusions and Outlooks - In this paper, we showed 
that both the event-by-event fluctuations of the IC and 
the continuous particle emission instead of sudden freeze- 
out largely modify the HBT correlation of produced pi- 
ons, so they should be included in more precise analyses 
of data. The IC fluctuations give smaller radii, without 
changing the dependence, which is a natural conse- 
quence of the presence of high-density spots at the early 
times. Continuous particle emission, on the other hand, 
does not change Rl but enhances the mr dependence of 
R s and inverts the mr behavior of R , which now de- 
creases with mr in accord with data. This is because, 
in this description, large-fcr particles appear mostly at 
the early stage of the expansion from a thin hot shell of 
the matter, whereas small-Zcy particles appear all over 
the expansion, and from larger portion of the fluid. The 
combination of these two effects can give account of the 
rriT dependence of Rl and R s and improves considerably 
the one for R a with respect to the usual version. 

We shall emphasize that these conclusions could 
only be reached because we have explicitly solved 3+1 
dimensional hydrodynamic equations and not simply 
parametrized the final flow as often done [2^. The re- 
sults are preliminary. In applying the CEM we had to 
make a drastic approximation, expressed by ea. l|ll|) . in 
order to make it feasible. It is likely that this is the 
reason why the discrepancy in R still persists. In addi- 
tion, there exist certainly dissipation effects. Since SPH 



is an effective description in terms of parameters i>i ap- 
pearing in eq.Jp), some smoothing of short- wave- length 
Fourier components is taken into account through the 
kernel W(x — Xi(t);h). We preferred not to explicitly 
include the viscosity at this stage, as it is still an open 
problem of hydrodynamics (see more details in Ref. 26]). 
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